Diamagnetic nature of stratified thin metals in visible range 
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It is numerically demonstrated that effectively strong diamagnetic resonance emerges at visible 
frequencies in stratified metal-dielectric metamaterials. The effective optical constants are extracted 
by two-complex reflectivity method. It is clarified that the effective diamagnetic response originates 
from local diamagnetism at stratified thin metals. The effective diamagnetism is crucially sensitive 
to the sturucture of unitcell. The effective diamagnetic response is always associated with effective 
plasma frequency and is therefore regarded as a magnetic component of the collective excitation. 
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I. INTRODUCTION 

Photonic metamaterials are trials to invent and pro- 
duce novel electromagnetic (EM) states by modifying pe- 
riodic structures in which constituents have the same 
optical constants with the bulk materials. The trials 
started at gigahertz, and has proceeded to terahertz and 
optical frequencies. The rapid progress was reviewed in 
[l|, l2J ■ The strategy is significantly different from nano- 
scale materials science in which quantum size effect on 
electronic states is a key. In metamaterials, novel effec- 
tive EM states are explored by controlling the geometri- 
cal structures. The effective EM states are solutions of 
Maxwell equations for a homogenized system obtained 
by coarse-graining and represent averaged EM fields ap- 
proximately. 

Recent metamaterial optics usually employs a retrieval 
way to evaluate effective optical constants Q- The pro- 
cedure is briefly explained as follows: When one knows 
complex transmissivity and reflectivity (called S param- 
eters) for a finitely thick object, one can numerically cal- 
culate (or retrieves) the effective permittivity and perme- 
ability by substituting the S parameters into the trans- 
mission and reflection formulas, which hold for a finitely 
thick material of effective optical constants. The way was 
first introduced in 2002 Q and has been called retrieval 
way. 

The retrieval way is algorithmic and not rigorous com- 
pared with the effective media description, which has 
been studied over one century 0, [H, @ since Maxwell- 
Garnet 0, H[ • The effective media descriptions are usu- 
ally derived for specified structures by assuming analyt- 
ical models and the long wavelength limit. The differ- 
ence between the effective media description and recent 
retrieval way was already described in the introduction 
of @. In short, the effective media description is derived 
analytically for ideally large objects under the long wave- 
length approximation, whereas the retrieval way trans- 
forms the S parameters into effective optical constants 
in a practical way. 
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In the retrieval way, there are two distinct points from 
the effective media description: (i) the interface of meta- 
materials play a crucial role because the S parameters 
are sensitive to the surface layers; (ii) the description us- 
ing effective optical constants turned out to be approx- 
imately valid for several artificial metamaterials beyond 
the long wavelength limit [9( . By use of the retrieval way, 
effective magnetic responses and negative refraction at 
optical frequencies have been demonstrated experimen- 
tally [13, El, G3 ■ Those metamaterials exhibit novel op- 
tical states, which are exceptions to the statement that 
permeability is unity at optical frequencies fl3j . 

A simple theoretical analysis presented a possibility of 
super-resolution by a thin metal film [14j | and the sub- 
diffraction-limited image was shown experimentally [l5l |. 
It is expected that near-EM-field enhancement is induced 
to compensate the loss at resonance with exploiting thin 
metallic layers. Another experimental result on transmis- 
sion was reported which supports the enhancement of EM 
fields by stratified thin metals • It was also shown that 
the enhancement is associated with high magnetic field 
in thin metals [T^. Besides, it was recently argued that 
negative refractive response appears in stratified stacks 
composed of metal-dielectric several layers [l8|,[lj|- Thus, 
the magnetic response in metallic thin films attracts great 
interest, connecting to novel optical phenamena. 

In this paper, stratified metal-dielectric metamaterials 
(SMDM) are specified to clarify the effective magnetic re- 
sponse at visible frequencies. Typical optical properties 
of SMDM have been shown in an earlier publication [2(| • 
Here, the effective magnetic response is focused on and 
explored with modifying the structure of unitcell. Figure 
Q] shows the schematic drawing of SMDM and the coor- 
dinate configuration. The metal is set to be Ag and the 
dielectric is MgF2 . The thickness of metal and dielectric 
are 15 and 60 nm, respectively. The periodicity is then 
75 nm. The thickness of Ag layers was set to ensure using 
dielectric constant of bulk silver in literature [2l[ . 

In Sec. UH two-complex reflectivity method (TCRM) 
are described. Numerical results are shown in Sec. IIII1 
and the implication is explored to elucidate the origin of 
effective magnetic response. More discussion and com- 
parison with other types of magnetism reported so far 
are given in Sec. IIVI 
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FIG. 1: Schematic drawing for structure of SMDM and co- 
ordinate configuration. Gray indicates metal (Ag) and white 
dielectric (MgF2). Blue bold line indicates the unitcell of 
SMDM. Incident light sheds on xy plane with s polarization 
(i.e. Ei n ||y) and incident angle 8. The thickness of surface 
layer is denoted by d. In this configuration, the two products 
of ey/fi x and fi x Hz are determined from two incident angles. 
Details are given in Sec. [TTJ 



II. TWO-COMPLEX REFLECTIVITY METHOD 

In this study, bulk SMDM is analyzed which is thick 
enough to eliminate transmission, and full components 
of e and /it tensors are extracted by TCRM. Although 
the scheme of TCRM was reported already [2(| , the full 
description is provided in the following. 

Since the structure of SMDM is uniaxial, effective ten- 
sors of permittivity e and permeability [i are assumed to 
be diagonal and uniaxial: 



(1) 



The optical axis is set to be the z axis, and then the two 
relations of e x = e y and \i x = /i y hold in Eq. ©. 

In linear and axial media, electric and magnetic flux 
densities are expressed using permittivity and permeabil- 
ity tensors as 
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where £o and Ho are permittivity and permeability in 
vacuum, respectively. 

First, the boundary conditions of Maxwell equations 
are taken into account. In Fig. [TJ incident wavenumvcr 
vector k; n is in the xz plane and then refracted wave has 
the wavevector k in the xz plane. We here introduce 
normalized refracted wavevector k = k/fco where ko is 
wavenumber of light in vacuum, and can write k as 



k = 



(4) 



Then, complex reflectivity r s under s polarization (that 
is, Ei n ||y) in Fig. Q] is easily derived from the Maxwell 
boundary conditions [22T ]: the equation for r s is some- 
times called Fresnel formula [23| . After simple modifica- 
tion, the next equation is obtained: 



k z (9) _ n cosdl -r s {9) 
H x no l + r s (6)' 



(5) 



where 9 is incident angle, and no denotes the refractive 
index in vacuum. In particular, k z (0) is refractive index 
along the z axis. 

Next, the equation of dispersion is included in TCRM. 
When the EM fields are monochromatic plane waves and 
proportional to exp(ik ■ r — iivt), the two of the Maxwell 
equations are 



ikxE = iwB, 
ikxH = — iwD. 



(6) 
(7) 



Note that the fields E, H, D, and B in Eqs. © and © 
are independent of (r, t). 

When the solution of plane wave is assumed, it is gen- 
erally allowed that the e and fi tensors depend on the 
frequency to and wavevector k; that is, the e and \i ten- 
sors have frequency and spatial dispersions (e = £(w,k) 
and fx = fj,(w, k)) [13| . However, it is for the present as- 
sumed that the e and [i tensors are independent of k. In 
solid crystals, the /z tensor is set to be unity at optical 
frequencies and independent of k, and the spatial dis- 
persion of e tensor is often so weak that it is frequently 
negligible [l3j . 

Equations © and J7J) yield Eq. ©, and Eq. © is 
derived with help of Eqs. © and ©: 
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Equation 
relation 



10]) is obtained by using k = fcpk and the 
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The determinant of coefficient matrix for E in Eq. (jXT 
has to be zero to have nontrivial solution; therefore, Eq. 
(fTU)l yields a general constraint between e and \x for a 
given k. The constraint condition is equation of dis- 
persion and can be regarded as generalization of Fres- 
nel equation. Historically, Fresnel equation was derived 
under the condition of \x = 1, gave the equation concern- 
ing e for a given k, and has been usually employed to 
explain optical phenomena, e.g. birefringence, in solid 
crystals |13j . 
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Under s polarization, Eq. (fT0|) is modified using Eq. 
© and Q : 
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We thus obtain the equation of dispersion under s polar- 
ization, 



k z (9) 2 , k x (0? 
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where k x (9) = no sin 6* in Fig. Q] 

Combining Eqs. ([5|) and (fTT|) . we obtain the next equa- 
tion: 



/ no cos 9 1 — r s (9) \ 2 n 2 , sin 2 i 
V /io l + r s (0)/ /i^Mz 



(12) 



When one knows the value of complex reflectivity r s by 
optical measurement or computation, Eq. (| 1 2[) is equa- 
tion for e y /fi x and fi x fi z . Two different angles (#i ^ 82) 
yield sin^i ^ sin #2, and enable to evaluate s y /fi x and 
fi x fi z uniquely. 

By permutating the configuration (x,y,z) — » (y,z,x), 
we can figure e z /fi v and fi y fi x - Then, incident light sheds 
on the yz plane and has the polarization of E in ||z. The 
plane of incidence is the xy plane. In this second step, 
Eq. (fT2|) is modified as 



n o cos0 1 - f s ((t))\ 2 | ngsin 2 ^ 
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FIG. 2: Effective optical constants of SMDM, extracted with 
the two incident angles of 0° and 15°. The surface thickness 
d of SMDM is 20 nm. (a) : Effective permittivity. Upper solid 
line denotes Ke(e x ) and upper dotted line does lm(e x )- Lower 
two lines represent e z likewise, (b): Effective permeability. 
Lines indicate fj, x and fi z with the notation similar to (a), 
(c): Effective refractive index along the x and z axes. The 
real parts are shown with solid lines and the imaginary parts 
with dotted lines. 
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where r s is complex reflectivity on the yx plane and <f> is 
incident angle in the xy plane. In evaluating s z /fx y and 
fiy fi x similarly to e y j n x and fi x fx z , it is required that the 
object should be semi-infinitely thick along the x axis. 

Since the relations of fi x = n y are assumed for SMDM, 
fi x is evaluated in the second step; succeedingly, the rest 
components of e x , e z , and fi z are uniquely obtained ex- 
cept for the sign. The sign of fi x is set to be fi x « 1 at 
800 nm because the fi x spectrum keeps constant and does 
not show any resonant behavior around the wavelength. 
The sign of fi x at different wavelength is taken to avoid 
discontinuous jump as shown in Fig. [5J 

Complex reflectivity r s was computed exactly, whereas 
f s was evaluated quite precisely, that is, within a-few- 
percent numerical error. The numerical way to compute 
r s was reported earlier (20| . The precision in the compu- 
tation for complex reflectivities ensures that the full com- 
ponents of tensors are determined within a-few-percent 
numerical error. 



III. NUMERICAL RESULTS 

A. TCRM Analysis 

Figure [5] presents the full components of (a) e and (b) 
fi tensors of SMDM of surface thickness d — 20 nm. The 
tensors are extracted by TCRM with the two angles of 0° 
and 15°. Other sets of two incident angles result in the 
same spectra of e x , s z , and n x within the numerical error. 
This consistency indicates that e x , e z , and fi x are inde- 
pendent of incident angle and consequently wavevector 
k. 

In TCRM analysis for SMDM using the angle pairs 
of (O°,0) and (0° ,</>), it is easily verified from Eq. 1(12]) 
that n z depends only on the 6(^= 0°). The three com- 
ponents of e x , e z , and n x could in principal depend on 
the incident angles; however, they have been found to 
be independent of incident angles and therefore to be 
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FIG. 3: Electric field vs the averaged field. Electric field is 
calculated at A p , e fl under normal incidence on the xy plane 
with s polarization. Solid line indicates forward propagation 
only, that is, incident and refractive components. Dashed 
line represents the electric field calculated by the effective 
refractive index n z in Fig. [2jc). Vertical dotted line indicates 
the interface. 
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well-defined for SMDM. On the other hand, the compo- 
nent fi z (9) is exceptionally evaluated for each non-zero 9. 
Thus, all the components of e and \i tensors are deter- 
mined self-consistently and in a well-defined manner at 
visible frequencies in Fig. [2] 

The e x component has typical dispersion of Drudc 
metal (upper solid and dotted lines in Fig. [Ha): the 
real and imaginary parts of s x , respectively), and the 
effective plasma frequency w PiC ff, which is defined by 
Re[e x (o;p ! cff)] = 0, is located at A PiG ff = c/27rw Pje ff = 468 
nm. On the other hand, the e z spectrum exhibits typi- 
cal dielectric properties for the incident light of E in ||/- on 
the xz or yz planes (lower solid and dotted lines in Fig. 
^ a): the real and imaginary parts of s z , respectively). 
Concerning permeability components, \i x is close to 1 at 
360-800 nm (upper solid and dotted lines in Fig. EJb): 
the real and imaginary parts of fi x , respectively), while 
\x z goes down to 0.57 at 477 nm and shows strong dia- 
magnetism (lower solid and dotted lines in Fig. [2jb) : the 
real and imaginary parts of respectively). Refrac- 
tive indexes along the x and z axes are displayed in Fig. 
[21c) ; the upper solid and dotted lines respectively denote 
Tie(n x ) and Im(n x ), and the lower solid and dotted lines 
respectively stand for Ke(n z ) and Im(n z ). The imaginary 
parts of n x and n z are small below A PjC ff , and correspond 
to efficient EM-wave propagation in SMDM. 

In Fig. [3l we test the validity of description by effective 
refractive index at A PjC ff. The interface of vacuum and 
SMDM is set at z = and incident light illuminates on 
the xy plane normally. Solid line denotes electric field of 
incident light in vacuum (z < 0); it represents refractive 
component of electric field in SMDM (z > 0). Dashed 
line presents the effective field profile calculated by using 
n z in Fig.^c) and reproduces the averaged field profile in 
SMDM approximately. Thus, TCRM extracts averaged 
EM response in SMDM. As was already confirmed in [20j |. 
the effective description has relevant physical meanings 
and works fairly well in the visible range of interest. 



FIG. 4: Comparison of (a) B z and (b) H z at 477 nm, the mini- 
mum of Re(ti z ). The region of vacuum (z < 0) shows incident 
light only, for simplicity; the incident angle is 6 = 15°. The 
region of SMDM (z > 0) contains both forwardly and back- 
wardly propagating components. Arrows indicate the vectors 
at typical points (dots). In particular, the directions of B z 
and H z in metallic layer (20 < z < 35 nm) are clearly oppo- 
site. 



B. Origin of effective diamagnetism 

As shown in Fig. [21 effective diamagnetic response ap- 
pears in the \i z spectrum. This implies that magnetic 
induction B z and magnetic field H z behave in a different 
way. Figure S] displays the actual distribution of (a) B z 
and (b) H z as image plot. Arrows denote the vectors 
at dotted points. The wavelength is at 477 nm, where 
Re(/i z ) is minimun. In vacuum (z < 0), the incident light 
travels with incident angle 15°; only incident component 
is shown for simplicity in Fig.[4j In SMDM [z > 0), both 
forward and backward components are displayed simul- 
taneously because reflectance R is not small (R ~ 0.7) 
and the backward components connected to inner multi- 
reflection is not negligible. The computation for B z at 
each point is carried out from Eq. ([6|), and on the other 
hand magnetic field H z is calculated from Eq. ([7|). 

Figure [4] obviously shows that the directions of B z and 
H z are opposite in metallic layer (20 < z < 35 nm). Di- 
electric layers behave like paramagnetic materials. Be- 
sides, the H z component in the metallic layer is far more 
intense than the H z in vacuum and dielectric layers. In 
this way, the thin metallic layers are found to be diamag- 
netic along the z axis. The diamagnetic nature does not 
come from electronic spin states of silver but from the 
superimposed EM fields peculiar to mesoscopic structure 
of SMDM. The backward component is indeed responsi- 
ble for the local diamagnetism. Thus, the effective dia- 
magnetic fj, z response is found to originate from the thin 
metallic layers. 

The diamagnetic response of metallic layer in SMDM 
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FIG. 5: Spatially dispersive behavior of \x z resonance. All the 
spectra are obtained by TCRM with the two angles of 0° and 
0°). However, as described in the text, \i z depends only 
on the in Fig.[T] The spectrum at 15° is the same with that 
in Fig. [U;b). The spectra at 9 = 30°, 45°, and 60° are offset 
for visibility by 0.5, 1.0, and 1.5, respectively. Horizontal 
dotted lines indicate Re(/i z ) — 1.0 for each spectrum. 



is confirmed numerically in Fig. |4j It is also easily derived 
from Maxwell equation. In Fig. |4j Eqs. (J6j> and Q in a 
metallic layer are explicitly written as 

iluB z = ik k x E y , (14) 
sgn(k z )ik k z H x - ik k x H z = -iuj£ Ag E y , (15) 

where sgn(fc 2 ) takes 1 or — 1, and denotes the forward and 
backward propagation for the z axis. From the computa- 
tion for actual EM distribution, \H Z \ is found handreds 
times larger than \H X \; i.e. \H Z \ 3> \H X \. Consequently, 
the first term in the left hand side of Eq. Q15p is negli- 
gible to a good approximation. Then, deviding Eq. (|14p 
by Eq. (fl~5|) . the sign of the real part of B z /H z turns out 
to be determined by a factor £A g : 

Rc(B z /H z ) oc Re(e Ag ). (16) 

At the wavelength in Fig.|H Re(eA g ) < directly results 
in Re(B z / H z ) < in the metallic layer. Thus, the dia- 
magnetic nature of thin metallic layers is confirmed from 
the simple calculation. 



C. Properties of \i z resonance 

Figure[5]shows the dependence of \x z on incident angle. 
The \i z spectra in Fig. [5] vary for incident angles and in- 
dicates spatial dispersion. Note that fj, z at off-resonance 
above 540 nm is almost independent of 9. The spatial 
dispersion is a resonant behavior. 

Figure [6] displays the effect of structure of unitcell on 
the effective \i z . Figures SHU are the results regarding the 
surface thickness d = 20 nm; the spectra is shown again 
in Fig. [51 From bottom to top, the thickness d varies from 
to 30 nm, respectively. By changing the surface thick- 
ness, the structure of unitcell in SMDM varys; see Fig.[T] 
SMDM with d = and 10 nm show diamagnetic prop- 
erty in a wide range. The increase of surface thickness 



d gradually reduces diamagnetism and finally eliminates 
the diamagnetism at d = 30 nm. Since the structure of 
unitcell directly affects the phase of EM wave in SMDM, 
it is crucially significant how to set the unitcell in re- 
alizing the effective diamagnetism. SMDMs of different 
unitcells are different metamaterials. Figure [6] is the ex- 
plicit indication. It is also consistent with the conclusion 
in a recent numerical study [24j , showing that the inter- 
faces play a vital role in photonic slabs of finite thickness, 
in extracting effective constants by the retrieval way. 

As is shown in Figs. [2] and O the effective diamag- 
netic response emerges at about A p>e ff, comes from the 
stratified metallic layers in SMDM. In order to observe 
the diamagnetic resonance in SMDM, it is essential how 
to select the structure of unitcell. As the A P]e ff shifts 
with varying the ratio of metal in the unitcell, the dia- 
magnetic resonance follows the shift of \ Pye g . It strongly 
suggests that the diamagnetic resonance is magnetic part 
of the effective plasma resonance. The spatially disper- 
sive behavior in Fig. [5] also shows a typical characteristic 
of resonance. 



IV. DISCUSSION 

From computational and theoretical studies, it was re- 
ported that effective magnetism takes place in the sys- 
tems composed of dielectric @, [25[ and metallic rods 
prj \2l\ , results from the Mie resonances peculiar to the 
rods, and is spatially dispersive [f||2t|. These magnetic 
responses stem from the local resonance due to geomet- 
rical shape of constituents. 

The effective diamagnetism in SMDM also inherited 
from the the local diamagnetic nature of thin metals, 
which was shown in Sec. IIII B[ There is a remarkable 
characteristic to the diamagnetism of SMDM, which is 
absent in the systems of the local Mie resonance. In the 
Mie systems it is known that the resonances always in- 
duce not a small amount of absorption loss, whereas the 
diamagnetic resonance is almost free from loss as seen 
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FIG. 6: Permeability spectra of SMDM with various unit- 
cells: the thickness d of surface layers from to 30 nm. All 
the spectra are computed at the incident angle of 15°, and 
are presented without vertical offset. The spectrum of d — 20 
nm in Fig. [5] is displayed again for comparison. 
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in Fig. Htc). Although small loss [i.e. small Im(n z )] ex- 
ists, it appears irrespective of the diamagnetism; indeed, 
the small loss also exists in the non-magnetic SMDM of 
surface thickness d = 30 nm in Fig. [5] [2(|. This kind 
of magnetic resonance with low loss has not been found 
so far to our knowledge and is a new class of magnetism 
in metamaterials. Consequently, it is expected that low- 
loss, effective magnetic resonance at visible frequencies 
will be realized based on this magnetism. 

Thick enough SMDM have been focused on in this pa- 
per. Of course, this system is not easily fabricated in 
reality. However, it is a good touchstone to examine how 
the effective optical responses emerge in periodic metal- 
dielectric systems. In the present analysis, there is no 
evidence to suggest negative refraction in SMDM, which 
was lately argued in finite systems [H, [l9[ . If negative 
refraction takes place in finitely thick SMDM, it would 
be due to an interference effect in finite systems. 

V. CONCLUSIONS 

The full components of effective tensors of e and fi in 
SMDM has been numerically extracted by TCRM in a 



well-defined way. As a main result, the effective strong 
diamagnetic resonance along the optical axis has been 
clarified. The resonance is spatially dispersive. It has 
been revealed that stratified thin metals in SMDM are 
the source of effective diamagnetism. The diamagnetism 
of thin metals results not from local Mie resonances but 
from the EM states realized in SMDM. The diamagnetic 
resonance is induced always near the effective plasma fre- 
quency; it can therefore be regarded as a magnetic com- 
ponent of the collective excitation. The diamagnetic na- 
ture of thin metals is experimentally feasible, and it will 
be a good test to measure nonlinear optical responses via 
diamagnetic resonance. 
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